Spitzer Transit and Secondary Eclipse Photometry of G J 436b 
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ABSTRACT 



We report the results of infrared (8 fim) transit and secondary eclipse pho- 
tometry of the hot Neptune exoplanet, G J 436b using Spitzer. The nearly 
photon-limited precision of these data allow us to measure an improved ra- 
dius for the planet, and to detect the secondary echpse. The transit (cen- 
tered at HJD = 2454280.78149 ± 0.00016) shows the flat-bottomed shape typ- 
ical of infrared transits, and it precisely defines the planet-to-star radius ratio 
(0.0839 ± 0.0005), independent of the stellar properties. However, we obtain 
the planetary radius, as well as the stellar mass and radius, by fitting to the 
transit curve simultaneously with an empirical mass-radius relation for M-dwarfs 
(M = i?). We find i?^ = = 0.47 ±0.02 in solar units, and Rp = 27, 600 ±1170 
km (4.33 ± 0.18 Earth radii). This radius significantly exceeds the radius of 
a naked ocean planet, and requires a gaseous hydrogen-helium envelope. The 
secondary eclipse occurs at phase 0.587 ± 0.005, proving a significant orbital ec- 
centricity (e = 0.150 ± 0.012). The amplitude of the echpse (5.7 ± 0.8 x 10"^) 
indicates a brightness temperature for the planet of T = 712 ± 36K. If this is 
indicative of the planet's physical temperature, it suggests the occurrence of tidal 
heating in the planet. An uncharacterized second planet likely provides ongoing 
gravitational perturbations that maintain GJ436b's orbit eccentricity over long 
time scales. 

Subject headings: Planetary systems - stars: individual (GJ436) - stars: low mass - 
stars: fundamental parameters - infrared: stars - eclipses 
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Introduction 



The transit and secondary eclipse of an extrasolar planet allow us t o deduce its physical 



prop erties to a de gree that is not po ssible in other observing geometries (jCharbonneau et al. 



20071 ). G J 436b flButler et al 



Neptune-sized planet ( iGillon et al. 



2004 



was recently discovered to be the first transiting 



2007al ). opening new parameter s pace for exop 



studies. It orbits an M-dwarf star lying 10 pc from our solar system (IManess et al. 



anet 



20071 ). 



In order to constrain the bulk composi tion and internal structure of transiting planets 



( ISeager et al. 



2007 



Fortney et al. 



20071 ). precise radii and temperature measurements are 
needed. The relatively small size of GJ 436A (~ 0.4 solar radii) enhances the planet-to-star 
contrast during transit and eclipse. Nevertheless, the shallow depth of the G J 436b transit 
(0.006) is a challenge for ground-based photometry. Altho ugh ground-based observers are 



achieving impressive levels of precision (jWinn et al. 



20071 ). photometry from space-borne 



platforms remains the preferred observational technique for the highest-precision transit 



have 


)een dominant 


(Charbonneau et al. 


2005; 


Demine et al. 


2005, 


2006; 


Harrington et al. 


2007; 


Knutson et al. 


2007 


)• 



In this Letter, we report Spitzer 8 fim transit and eclipse observations of G J 436b, 
and we use these data to refine estimates of the planet's radius, temperature, and internal 
structure. 



2. Observations 



The announcement of G J 436b transits (IGillon et al.ll2007al ) was fortuitously concurrent 
with a window of observability using Spitzer. Accordingly, we immediately scheduled 
observations of one transit, and one secondary eclipse, under our GO-3 Target of 
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Opportunity (ToO) Program (J. Harrington, P.I.). Since the precision required for 
measurements of this typ e is daunting, observati ons must be carefully designed to limit 



instrumental svste matics (IHarrington et al 



the G J 436b orbit (IManess et al. 



20071 ). Moreover, the reported eccentricity of 



20071 ) adds significant uncertainty to the timing of the 



secondary eclipse observations. Our community ToO program thus works with cooperating 
teams to design observations and analyze data in line with the best practices gleaned from 
experience. We solict collaborations from cooperating teams that discover suitable targets. 



Both observational sequences for GJ436 used the IRAC instrument (iFazio et al 



200J) 



in subarray mode, at 8 fim only. The transit sequence consisted of 0.4-second exposures in 
blocks of 64, obtaining 445 blocks (204 minutes). The secondary eclipse sequence was the 
same, but used 780 blocks (356 minutes). We planned the eclip se observations based on 



10^ bootstrap trial fits to the Doppler data (IManess et al. 



20071 ). to define the probability 



distribution of eclipse time. Because of the well kno\ yn ramp-up in the sensitivity of the 



IRAC 8 /im detector during long observing sequences (iKnutson et al. 



2007 



Harrington et al. 



20071 ). we offset the transit observations to begin ~ 2 hours before transit center. 



3. Data Analysis 



3.1. Photometry 



Because GJ436 is bright at 8 /im, and the zodiacal background is weak in comparison, 
simple aperture photometery attains nearly photon-limited precision. Our photometry first 
applies the calibration information contained in the FITS headers, to convert the signal 



levels to electrons. Within each 64-frame block, we drop t 



le first frame and the 58th 



frame, due to known instrumental effects ( IHarrington et al 



2007 



Knutson et al 



2003). We 



examine the time variation of signal level for each pixel in the remaining 62 frames, and 
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correct pixels in frames that are discrepant by > 4cr to the median value for that pixel 
(this removes energetic particle hits). We sum the intensity in an 8- by 8-pixel square 
aperture centered on the star in each frame, including fractional pixels, and sum again 
over the 62 frames in the block. We varied the aperture size to verify that an 8-pixel box 
produced the lowest noise, but this dependence is not strong. We fit a Gaussian to the peak 
in a histogram of pixel intensities for each block to determine, and subtract, the average 
background level. We used the same background value for all 62 frames in a block. 

We calculated the expected noise level for the photometry, based on the Poisson 
electron counting noise (dominant), and read noise (small). Comparing the aperture 
photometry for the 62 frames within each block, we find that these photometry errors are 
distributed as Gaussian noise, with a dispersion merely 3.5% greater than predicted. The 
intensities for the 445 transit blocks are illustrated in Figure 1, top panel. We examined 
the block-to-block variation in intensity for these points after removing the best transit 
fit, and we find a Gaussian distribution, with a standard deviation of 7 x 10~^. Since we 
detect ~ 3.1 X 10^ electrons per block (Figure 1), we expect a photon-limited precision of 



5.7 X 10 We thus attain about 80% of photon-limited S/N, consistent wit 



Spitzer photometry at this wavelength (IKnutson et al. 



2007 



Harrington et al. 



1 pre vious 



20071). 



IRAC photometry at 8 fim is known to exhibit a gradual 



in sensitivity, due to filling of charge traps in the detectors (IKnutson et al 



Harrington et al 



ly increasing ramp- 



up 



20071 : 



20071 ). This ramp is visible in the top panel of Figure 1, but is weaker than 



usual for the transit data (the ramp varies due to prior usage of the detector). We removed 
it by masking out the data near transit, and fitting a parabola to the out-of-transit points. 
We have considerable experience in fitting to this ramp, via our monitoring of GJ876 
(program 30498). Even strong ramps can be fit by the sum of a linear plus logarithmic 
function, using linear regression. We applied this more elaborate procedure to the GJ436 
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transit ramp, but found no significant difference with the simple parabola fit. In the case of 
the secondary eclipse (Figure 2), the ramp is stronger, but is still well removed by our full 
linear + logarithmic fit. We conclude that this ramp is properly reproduced in both cases, 
and does not contribute significantly to our errors. 



3.2. Transit Parameters 



lar limb darkening. Not 



Richardson et al. 



fl2006h 



A feature of IR transit measurements is the virtual lack of ste^ 
only does this produce a simple box-like shape for the transit, but 
suggest that it can increase the radius precision for a given level of photometric precision. 
Our analysis adopts the (small) limb darkening for GJ436A based on a Kurucz model 



atmosphere for 3500/5.0/0.0 in Teff/log(g)/[M/H]. We verified that changing the ste . 



ging the ste 


lar 


IBean et al. 


2006) 



temperature, gravity, or metallicity within the errors (iManess et al 
has negligible effect, because the limb darkening remains small over the plausible range. We 
integrated the stellar center-to-limb intensities in the Kurucz model over the bandpass of 
the IRAC 8 fim filter to obtain the limb darkenin g appro p riate to this IR transit. Since 
this small IR limb darkening is not included in the I Claret I (l2000l ) prescriptions, we generate 
theoretical transit curves numerically. 

We compute theoretical transit curves by tiling the star in a latitude-longitude grid 
with zone spacing of 0.18 degrees, and applying the IRAC 8 fim limb darkening. We pass 
the planet across the numerical star in steps of 0.01 stellar radii, with the planet radius and 
impact parameter specified in units of the stellar radius. To increase precision, stellar zones 
at the edge of the planet are adaptively sub-sampled in a 10 x 10 finer grid . We v erified the 



Mandel and Agol 



(120021 1 analytic 



code's precision (better than 10~^) by comparing to the 
non-linear limb darkening cases, and by comparing the depth of synthetic transits to the 
planet-to-star area ratio {Rp^/R^'^), for the case when limb darkening is identically zero. 
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Fitting to high 



precision transit p h otometry requires a d etermination, or assumption 



of the stellar mass (iBrown et al.ll200ll ). 
for GJ436A, based on th e obse r ved 
mass-radius relation from 



Gillon et al. 



(l2007a| ) adopted 0.44 solar masses 



uminosity, and they cited the empirical M-dwarf 
RibasI (j2006l ) {R = M in solar units), to justify 0.44 solar radii for 
the stellar size. Our fit procedure is somewhat different. Given the lack of limb darkening, 
we can immediately determine the ratio of planet to stellar radius as 0.0839 ± 0.0005 from 
the depth of the transit (Figure 1). With this value fixed, we generate a grid of transit 
curves for a range of impact parameters. At each impact parameter, we vary the adopted 
stellar mass, and compute the transverse velocity of the planet acros s the star. This 



computation uses the orbital elements from a fit to the Doppler data (iManess et al.l 120071 ). 
constrained by the secondary eclipse time (see below). We vary the stellar radius to convert 
the radius increments on the abscissa of the synthetic transit curve to orbit phase, using the 
calculated transverse velocity. We include a shift i n phase for the synthetic transit curve. 



to allow for imprecision in the 



Gillon et al. 



(l2007al ) ephemeris. In this manner, we find the 



best fit stellar radius versus stellar mass, and a revised transit time. We estimate the stellar 
radius precision from the variation in at a given mass. 

We intersect the radius versus mass relatio n (R ~ M°'^ ^) from the fitting procedure 



20061 ) to find the best stellar mass 



with the empirical mass-radius relation R = M ( iRibasI 
and radius, and planet radius, at each impact parameter. Repeating this over a grid of 
impact parameters, we adopt the best fit from the global minimum x^- We determine the 
error range from Sx^, and from visually inspecting the quality of the fits, paying particular 
attention to ingress/egress. Our fitting always uses the unbinned data (Figure 1, top), but 
we bin the data for the lower panel of Figure 1, to better illustrate the quality of the fit. 
The derived time of transit center is HJD = 2454280.78149 ± 0.00016. 
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3.3. Secondary Eclipse 

The secondary eclipse is shown in Figure 2. The top panel plots the bulk of the data 
(omitting some points at the outset); the eclipse occurs near the end of the observational 
sequence, at phase 0.587 ± 0.005, with amplitude 5.7 ± 0.8 x 10~^ in units of the stellar 
intensity. Like the transit, all of our fits to this event were made on the original, unbinned 
data. However, for clarity, the lower panel of Figure 2 shows binned data, expands the 
phase scale, and overplots the best fit eclipse curve. Our fit constrains the duration of 
eclipse to equal the duration of transit, finding only the amplitude and central phase. 

Table 1 summarizes our results for transit and secondary eclipse. 



4. Results and Discussion 



Our result for the stellar mass and radius is M = R = 0.47 ± 0.02 in solar units. 
We are encouraged that these are close to values (0.44) constrained by independent data 



(IManess et al. 



20071). Our derived planet radius is Rp = 27, 600±1170 km (4.33±0.18 Earth 



radii). We c onclude that this p lanet is larger than originally indicated by ground-based 



photometry (IGillon et al. 



exclusively solid composition (ISeager et al. 



2007a[). Since this radius is significantly lar ger than all planets of 



2007 



Fortney et al 



2003), G J 436b must have 



a significant, gaseous, hydrogen-helium envelope. After our transit analysis was complete, 
we became awa r e of ari independent analysis of the transit data (but not the eclipse data) 



by 



Gillon et al. 



(l2007bl ). These authors do not vary the stellar mass in their fit, but they 



obtain a very similar radius for this planet, and arrive at essentially the same conclusion. 

The observed phase of the secondary eclipse, (f> = 0.587 ± 0.005 indicates that the 
orbit of G J 436b is significantly eccentric. Assuming a longitude of pericenter m = 0, 
the magnitude of the observed timing offset indicates a minimum orbital eccentricity. 
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0.137 ±0.007. 



U sing; the constraints provided by the observed times of central transit, = 2454222.616 



HJD f Gillon et al. 



2007ah . our = 2454280.78149 HJD, and the observed secondary 



echpse at Tg = 2454282.33 ± 0.01 HJ D, we obtained a set of single-planet Keplerian fits to 



the radial velocity dat a published 



resampling procedure (jPress et al 



Maness et al. 



(120071 ). A straightforward bootstrap 



1992) yields e = 0.150 ± 0.012, = 343 ± 14°, and 



M = 0.070 ± 0.003 Mjup. 

As a consequence of its non-z ero orbital eccentricity, G J 436b is likely experiencing 



asynchronous rotation. 



Hut 



(Il98ll ) gives an expression for the spin period of a zero-obhquity 



spin pseudo-synchronized planet: 



spin 



:i + 3e^ + |e^)(l-e^)^/^ 
1 + f e2 + f e4 + 



orbit 



(1) 



For GJ436b, we find Pgpin = 2.32 d, which yields a 19-day synodic period for the star as 
viewed from a fixed longitude on the planet. The large orbital eccentricity also indicates 
that a significant amount of tidal heating must be occurring. To second order in eccentric ity, 



the tidal luminosity of a spin-synchronous planet (jPeale fc Cassen 
given by: 

dE _ 21 k2 GMlnRle^ 



1978 



Mardling 



20071) is 



dt 



2 Q 



(2) 



where k2 is the planetary potential Love number of degree 2, n is the orbital mean motion, a 
is the orbit al semimajor a x is, an d Q is the planet's effective tidal dissipation parameter. The 



analysis of 



Levrard et al. 



(120071 ) indicates that the tidal luminosity of an asynchronously 
rotating planet with e ~ 0.15 will exceed the value implied by the above expression by a 
small amount. 

If we adopt Teg = 33507^ for GJ436A, take a zero albedo for the planet, and assume 
a uniform re-radiation of heat from the entire planetary surface, we obtain a planetary 
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Teq — 642 K. The somewhat higher temperature (T = 712 ± 36K) imphed by the secondary 
echpse depth could arise from inefficient transport of heat to the night side of the planet, 
from a non-blackbody planetary emission spectrum, from tidal luminosity, or from a 
combination of the three. 

We can estimate Q by assuming ^4 = 0, uniform re-radiation, and a blackbody 
planetary SED to find a fiducial tidal luminosity of 4.7 x 10^^ ergssec"^ and £ = 2.1 x 10^. 
Assuming k2 — 0.34 (the Jovian value), gives Q — 7140. This value is in rough accord 
with the Q- values measured for Uranus and Neptune. Banfield and Murray (1992) derive 
1.2 X 10"^ < (5jv < 3.3 X 10^ for Neptune, whereas Tittemore and Wisdom (1989) employ 
the Uranian satellites to derive Qu < 3.9 x 10^. If the planet has maintained e ~ 0.15 for 
biUions of years, Q ~ 7000 indicates that the planet has radiated tidal energy comparable 
to the orbital energy and in excess of 100 times its own gravitational binding energy. 

Furthermore, Q ~ 7000 imphes a circularization timescale ~ 3.0 x 10^ yr. Indeed 
for any range oi Q - which is uncertain even for solar system bodies - the circularization 
timescale is < 10^ years. It is thus highly likely that an as-yet uncharacterized second 
planet is providing ongoing gravitational perturbations that allow GJ 436b's eccentricity to 
be maintained over long timescales. 

We are grateful to the staff at the Spitzer Science Center for their prompt and efficient 
scheduling of our observations. 

Facilities: Spitzer. 
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Fig. 1. — Spitzer photometry of the GJ 436b transit with fitted transit curves. Top: photom- 
etry before basehne correction. The dashed hne is the adopted basehne. Bottom: Basehne- 
corrected data, binned to approximately 2- minute time resolution (137 sec), with the best 
fit transit curve. Note the flat bottom that proves a non-grazing transit. 
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Fig. 2. — Secondary eclipse photometry of GJ436b. Top: photometry showing the ramp in 
intensity, with the echpse marked at phase 0.587. The dashed hne is the adopted basehne 
for points with phase > 0.52. Bottom: Binned data (102-sec time resolution) shown in 
comparison to the best fit secondary-echpse curve, whose amphtude is 5.7 ± 0.8 x 10~^. 
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Table 1. Derived parameters for GJ436. 



Parameter Value 



Stellar radius'' 
Stellar mass 
Planet radius 
Impact parameter 
Transit time 
Orbit semi-major axis 
Orbit eccentricity 

Secondary eclipse phase 
Secondary eclipse amplitude 
Planet brightness temperature 
Planet mass 



0.47 ±0.02 
0.47 ± 0.02 
27, 600 ± 1170 km 

n Qi;+0.03 
U.OO_() Q2 

HJD = 2454280.78149 ± 0.00016 

0.0291 ± 0.0004 AU 

0.150 ± 0.012 

13.2 ±0.6 

0.587 ±0.005 

5.7 ±0.8 X 10-4 

712±36K 

0.070 ± 0.003 Mjup 



"■Radius constrained to equal mass, in solar units. 



